The nondeterministic Nagel-Schreckenberg traffic model with 

open boundary conditions 



S. Cheybani a,fe and J.Kertesz 6 ' c and M. Schreckenberg a 
a Theoretische Physik, Gerhard- Mercator Universitat, D-^7048 Duisburg, Germany 
b Department of Theoretical Physics, Technical University of Budapest, H-llll Budapest, 

Hungary 

Laboratory of Computational Engineering, Helsinki University of Technology, FIN- 021 50 Espoo, 

Finland 



Abstract 

We study the phases of the Nagel-Schreckenberg traffic model with open 
boundary conditions as a function of the randomization probabilities p > 
and the maximum velocity v max > 1. Due to the existence of "buffer sites" 
which enhance the free flow region, the behaviour is much richer than that 
of the related, parallel updated asymmetric exclusion process (ASEP, v max = 
1). Such sites exist for v max > 3 and p < p c where the phase diagram is qual- 
itatively similar to the p = case: there is a free flow and a jamming phase 
separated by a line of first order transitions. For p > p c an additional maxi- 
mum current phase separated by second order transitions occurs like for the 
ASEP. The density profile decays in the maximum current phase algebraically 
with an exponent 7 ~ § for all v max > 2 indicating that these models belong 
to another universality class than the ASEP where 7 = | 



I. INTRODUCTION 



Asymmetric exclusion processes (ASEP) play an important role in non-equilibrium statistical 
mechanics. The one- dimensional ASEP is a lattice model which describes particles hopping 
in one direction with stochastic dynamics and hard core exclusion. It was first introduced 
in 1968 to provide a qualitative understanding of the kinetics of the protein synthesis on 
RNA templates [I]]. It turned out, however, that - despite its simplicity - there are numerous 
further applications of the ASEP on the field of interface growth, polymer dynamics, and 
traffic flow §- §. 

Unfortunately, as far as traffic is concerned, the ASEP yields rather unrealistic results, 
because essential phenomena like acceleration or slowing down cannot be reproduced in this 
model. As a consequence, Nagel and Schreckenberg developed an extension of the ASEP 
resulting in a one-dimensional probabilistic cellular automaton model ||. According to the 
Nagel-Schreckenberg model the road consists of a single lane which is divided into L cells 
of equal size numbered by i = 1, 2, . . ., L and the time is also discrete. Each site can be 
either empty or occupied by a car with integer velocity v = 0, 1, . . ., v ma x- All sites are 
simultaneously updated according to four successive steps: 

1. Acceleration: increase v by 1 if v < v max . 

2. Slowing down: decrease v to v = d if necessary (d: number of empty cells in front of the car). 

3. Randomization: decrease v by 1 with randomization probability p if p > 0. 

4. Movement: move car v sites forward. 

Either ring (periodic boundary conditions) or open (open boundary conditions) geometry 
is considered. In the case of ring geometry cars move on a ring and the car density in the 
system keeps constant. Open systems, on the other hand, are characterized by the injection 
(extinction) rate a (/?), that means, by the probability a ((3) that a car moves into (out of) 
the system. 

For the maximum velocity v ma x = 1 the model is identical with the ASEP with parallel 
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update 0- H which has been solved exactly with periodic boundary conditions and 
recently also with open boundary conditions |J, [Kj. In this special case three regimes (free 
flow, jamming, and maximum current) can be distinguished from each other. The transition 
from the free flow to the jamming phase at the a = f3 - line for a, ft < l-yfp is of first order. 
The transition from the free flow (jamming) to the maximum current phase is continuous 
and takes place at the injection (extinction) rate a c (/3 C ) with a c (p) = /3 c (p) = l-\z1>- 
As it is common for traffic simulations cars are updated in parallel in the Nagel- 
Schreckenberg model, too, because this update scheme is the only one which models the 



formation of spontaneous jams occurring in real traffic [plj| . Systems with parallel update 



are furthermore characterized by strong short-range correlations, and therefore, short-range 



correlation functions play an important role here ||12|| , [f24| . 

Most of the work dealing with the Nagel-Schreckenberg model for v max > 1 impose periodic 
boundary conditions [|13j- p3fl . Much attention has been paid to the question of the tran- 
sition from the free flow to the jamming regime. According to the state-of-the-art this is a 
crossover rather than a sharp transition ]TB[- |2"5fl . Systems with periodic boundary condi- 



tions are furthermore characterized by a trivial density profile p(i) = p with 1 < i < L due to 
translational invariance. In this context it should be mentioned that short-range correlation 



functions are well-suited for the description of the free flow - jamming transition |2J]: The 
free flow regime is characterized by anticorrelations around a propagating peak, that is, in 
free flow cars are surrounded by empty space. At the critical density p c the anticorrelations 
are maximally developed, and for higher densities they vanish. Simultaneously, a jamming 
peak develops according to the fact that the back car is strongly slowed down in a jam. 
In the following, systems for maximum velocities v max < 10 are investigated for the more 
realistic case of open boundaries. Boundary conditions are defined as in |[26|| : At site i = 
0, that means out of the system a vehicle with the probability a and with the velocity v = 
Vmax is created. This car immediately moves according to the Nagel-Schreckenberg rules. If 
the velocity of the injected car on i = is v = (because site i = 1 is occupied by another 
car or because the front car is on site i = 2 and the injected car is slowed down by 1 due to 



randomization) then the injected car is deleted. At i = L+l a "block" occurs with proba- 
bility 1-/3 and causes a slowing down of the cars at the end of the system. Otherwise, with 
probability (3, the cars simply move out of the system. In [^] systems with open boundaries 
have been already analyzed for the randomization probability p = 0, i.e., by ignoring the 
randomization step. The most interesting feature of the deterministic Nagel-Schreckenberg 
model for maximum velocities v max > 3 and open boundaries is the existence of so-called 
buffers: As a consequence of the parallel updating and the hindrance an injected car feels 
from the front car at the beginning of the system spaces larger than v max develop between 
two neighbouring cars for high injection rates. This can be easily demonstrated by consid- 
ering i = and the first sites of the system i = l, • • •, 5 for a = (3 = 1 and v max > 3. At an 
arbitrary time t = t a car with velocity v max is injected on i = 0, i.e. 

injection : v max \ . . 2 . . 

After application of the NaSch rules on the system we have injected on i = 0, i.e. 

movement : . | . 2 . . . 

Correspondingly, we get for t = t + 1 

injection : v max | . 2 . . . 
movement : . | 1 . . . 3 

and for t = to + 2 

injection : v ma;r | 1 ... 3 
movement : | . . 2 . . 

As the car on site i = cannot move it is deleted and the situation starts over in the next 
time step: 

injection : v max | . . 2 . . 
movement : . | . 2 . . . 
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and so on. Obviously, one car is lost out of three injection possibilities. In the system, far 
from the boundaries, the distance between two neighbouring cars turns out to be alternately 
di = Vmax and d.2 = 2{v max — 1). That means, in addition to the expected v max empty sites 
larger gaps occur in the a — > 1, /3 — > 1 limit. We call these additional sites "buffers" because 
they have a buffer effect at the end of the system: Due to these sites the development 
of jamming waves is suppressed even for (3 < 1. The transition from the free flow to the 
jamming phase is of first order and accompanied by the collapse of the buffers. The effect 
resulting from the buffers do not depend on the maximum velocity if v m ax > 5 (for v max 
= 3,4 the buffer effect is not so strong as the buffers are not completely developed for that 
case). 

However, randomization is indispensable for the analysis of real traffic as it takes human 
behaviour into account: The behaviour of a car driver is not like that of a machine but 
rather contains unpredictable elements. In traffic, over-reactions when slowing down can 
be found as well as delays when accelerating, furthermore fluctuations when following a car 
(follow-the-leader situation) and so on. 

Besides this motivation it is of interest to compare the general p > 0, v m(M . > 1 case with 
the previously investigated models (p > 0, v ma x = 107 13 an< ^ P = ^> Vi > 1 P^fl)- The 
presented results were obtained by simulating a L = 1024 sample with at least 1000 runs with 
10 4 time steps each. In order to investigate the influence of randomization on the system 



we proceed similarly to |26]: Section II considers the behaviour of current and average 
occupation number in the middle of the system. Section III deals with density profiles, 
Section IV with short-range correlation functions. Finally, the results are summarized in 
Section V. 



5 



II. CURRENT AND AVERAGE OCCUPATION NUMBER AT THE MIDDLE OF 

THE SYSTEM 



The phase diagram for systems with probability p = 0.5 and maximum velocity v max = 2 
and v m ax > 5 is shown in Figs la,b (the case v max = 3,4 is similar to the latter case). For 
v mai = 2 the phase diagram (Fig la) is qualitatively the same as for the case v max = 1 ||, 
||10|| : The free flow and the jamming regime are divided by a straight line and for a > 0.35 
and j3 > 0.8 the system is in the maximum current phase. Obviously, the maximum current 
regime is smaller than for v max = 1 and there is no symmetry along the a=/3-line. For v max 
= 3 the maximum current regime is even smaller ( for a > 0.35 and (5 > 0.85 ) and the free 
flow/jamming border shows a slight bending. These tendencies are even stronger developed 
for higher maximum velocities v max with a maximum current regime for for a > 0.35, f3 > 
0.87, and v max = 4 (for a > 0.35 and (3 > 0.89, and v max > 5, see Fig lb). 
Another interesting feature of the nondeterministic case is that the course of the free 
flow/jamming border is totally different from that for p = 0. In Fig lc it can be clearly seen 
that this difference is a consequence of the vanishing buffer effect due to randomization. For 
randomization probabilities p > p c (p c = 0.1172 ± 0.008 for v max = 5) there is no sign of 
the buffer effect any more, and a (rectangular) maximum current regime develops for a > 

C^FreeFlow &nd f3 > f3 Jamming- 

In order to understand the nature of the transition between the phases we consider the 
average occupation number on the site i = ^-,p(i = ^), as proposed in p. Fig 2 shows 
p (i = tj) for Vmax = 5 as a function of the injection and the extinction rates (the average 
occupation number on i = ^ for any v max > 1 behaves similarly). It turns out that phase 
transitions in systems with maximum velocity v max > 1 show the following features: The 
transition from free flow to jamming is of first order indicated by a jump in p (i = ^). At 
the maximum current /free flow and the maximum current /jamming transition there is a 
jump in the derivative of p (i = ^) what is a hint at a continuous phase transition. These 
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conclusions will be confirmed in the next section which deals with the investigation of the 
corresponding density profiles. 

In the following we analyze the current due to the influence of the 

- boundaries 

- maximum velocity v max 

- randomization probability p 

The best way to investigate the influence of the left boundary is to consider the case j3 = 
1 for p = 0.5 where cars simply move out of the system. In Fig 3a free flow and maximum 
current phase can be clearly distinguished from each other. For maximum velocities v max 
> 5 the curves are nearly the same with a maximum at a ~ 0.35 becoming stronger with 
increasing v max (the occurrence of the maximum will be explained below). The dependence 
of the current on a and [3 for v max = 2, 3, 4 is qualitatively the same as for v max = 1. 
For the investigation of the influence of the right boundary we consider the case a — 1 for p 
= 0.5 (Fig 3b). Here, the current for v max > 5 does not depend on the maximum velocity. 
Furthermore, it seems to increase monotonously with increasing (3 also in the maximum 
current phase {j3 > 0.89). Investigations for system sizes L > 4096, however, show that the 
latter observation is just a finite size effect and that the current for v max > 5 and (3 > 0.89 
is constant. Apart from this, the curves for the current do not change in an essential way 
with increasing system size L, and therefore it is sufficient to investigate systems with L = 
1024 in the following. 

From the observations so far we can conclude that - as for the deterministic case - the 
behaviour of the system only negligibly changes when maximum velocities v max > 5 are 
considered. Therefore we restrict ourselves to the case v max = 5 in the following observa- 
tions. 

We will now investigate the influence of the randomization probability on the behaviour 
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of the system. It can be seen from Figs 4a,b that the buffer effect observed for the deter- 
ministic case vanishes with increasing p: For (3 — 1 in Fig 4a the maximum at a ~ 0.81 
resulting from the existence of the buffers moves to the left and becomes weaker and weaker 
(In Fig 4a we make an exception and consider the current at v max = 10 instead of v max = 
5, because for the latter case this effect is nearly invisible). For the injection rate a = 1, on 
the other hand, the buffer effect vanishes as soon as randomization probabilities p > p c are 
considered. 

To sum it up it can be said that as a consequence of the buffer effect, the course of the 
current in the maximum current phase deviates from the expected (constant) behaviour for 
(3 — 1, p = 0.5, and maximum velocities v ma x > 5 showing a slight maximum at a ~ 0.35. 
Besides, there are strong indications that a continuous transition from the free flow (jam- 
ming) to a maximum current phase develops with increasing randomization probability on 
the (3 — 1 - line (a = 1 - line). More convincing arguments for the existence of a maximum 
current phase, however, will be given in the following sections. 

III. DENSITY PROFILES 

Our observations so far consider the behaviour of the whole system and of the site i = ^. 
For the analysis of what happens on the other sites it is useful to investigate the density 
profiles. Of special interest in this context is the question in how far the density profiles 
reflect the transition between the phases. For that purpose we consider the density profiles 
for (3 = 1 (— > transition from free flow to maximum current), for (3 = 0.7 (— > transition from 
free flow to jamming), and for a = 1 (— > transition from maximum current to jamming). It 
turns out that - as in the case of v max — 1 - the free flow (jamming) phase can be divided 
into the regime AI and All (BI and BII). The following investigations are confined to the 
randomization probability p = 0.5, L = 1024 and v max = 5 (for v ma x > 5 and L > 1024 the 
density profiles are qualitatively the same). 

In Fig 5a the transition from free flow to maximum current for (3 = 1 can be clearly seen. 
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The free flow regime is characterized by oscillations at the beginning of the system dying 
out for i > 100 (if p = 0.5 and v max — 5) due to randomization and the density profile 
becomes constant. It can therefore be said that randomization blurs the influence of the left 
boundary. In the maximum current regime we do not have any oscillations at all. Instead, 
an analytic decrease of the density is observed becoming stronger with increasing a. This 
phenomenon can be easily understood hinder each other at the beginning of the 

system for high injection rates: The higher the injection rate the stronger the hindrance. 
As a consequence the density profiles in the maximum current phase do not depend at all 
on a in the middle and at the end of the system. At the beginning of the system, however, 
the density profiles decay as i~ 7 with 7 « 0.66 what is valid for all v max > 1. We conjecture 
that 7=| because the exponent converges to this value with increasing system sizes. That 
means that the cases v max = 1 and v max > 1 belong to different universality classes since 
the corresponding density profiles for the ASEP with parallel update decay as i -7 with 7 = 
\ at the beginning of the system ||, Jit]]. 

The free flow/maximum current transition is nicely reflected by the density at the end of 
the system: p(i=L) is proportional to the injection rate when the cars move freely up to a c 
= 0.35 and becomes constant in the maximum current regime. 

The situation is different when density profiles at the free flow - jamming border for (3 = 
0.7 are considered (Fig 5b). It can be easily seen that the transition is of first order as 
the density profile at a c = 0.278 is linear which is a typical feature of a first-order phase 
transition (see |TlJ] and references therein). 

The density profiles for the injection rate a = 1 are shown in Fig 5c. The course of the curves 
for [3 < 0.89 (/3 > 0.89) is typical for the maximum current (BII jamming) phase with an 
algebraic (exponential) decay at the beginning of the system due to the hindrance already 
described for (3 = 1. The decay becomes weaker with decreasing (3 and finally vanishes 
as the repercussion resulting from the blockage at i = L+l increasingly superimposes the 
hindrance effect at the beginning of the system. For extinction rates 0.6 ^ ^ 0.9 there is 
a slight increase at the end of the system which indicates a hindrance due to the blockage. 
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It just remains a border effect, however, and is of no relevance for the considerations in this 
article. 



IV. CORRELATION FUNCTIONS 



In this chapter we consider short-range correlation functions 



C(i,t) = < 77(i',t') r/(i+i',t+t') > v v - < T/(i',t') > r 



where 



r/(i',t') = 1 if site i' is occupied at time t ' 
r/(i',t') = else 

< ... >j5 describes the spatial and temporal average over all L sites i' and over times t' 
taken from our simulation of the steady state. 

The correlation functions are measured in the middle of the system where the influence of 
the boundaries is minimal. We do not only investigate the cases (3 = 1 (influence of the left 
boundary, Fig 6a) and a = 1 (influence of the right boundary, Fig 6b), but also (3 = 1-a 
(Fig 6c). For the latter case there are similar conditions at both boundaries and the system 
can be compared at best with the corresponding system with periodic boundary conditions. 
Therefore it is no surprise that the correlation functions in Fig 6c are qualitatively the same 



as those for periodic boundary conditions [E| 



What is interesting, however, is that a classification into free flow, maximum current, and 
jamming cannot be done when short-range correlation functions are considered. Instead, 
due to Figs 6a-c three regimes can be distinguished from each other: 

(a) Free Flow: The free flow is characterized by anticorrelations around a 
propagating peak at i = v max (t-l) with a shoulder at i = v ma x t, that is, 

in free flow moving cars are surrounded by empty space 

(b) Coexistence Regime (JI + Maximum current): The coexistence of 
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free flow and jamming manifests itself in the double peak structure of the 
correlation function. The jamming causes a maximum at i = -1 according 
to the hindrance the back car feels in the jam. 
(c) Jamming II (JII, "Superjamming"): The propagating peak disappears 
as a consequence of the fact that in free flow moving cars do not exist any 
longer. 

As it is obvious from the previous sections the transition from JI to JII is not a phase 
transition and does not change the behaviour of the system in an essential way. 
In correspondence with |^4| the critical injection (extinction) rate a C 2 {(3 C 2) for the JI- JII 
transition is defined by the vanishing of the propagating peak and takes place at /3 C 2 ~ 0.65 
(see also Fig lb). Unfortunately, an exact value for the Free Flow- Jamming transition can 
neither be given. As a consequence of the randomized oscillations in the density at the 
beginning of the system it is not possible to determine the critical injection (extinction) rate 
at which the influence of the right boundary reaches the left boundary. This is a significant 
difference to the deterministic case [EBJ where the oscillations of the free flow phase form a 
well-defined pattern due to the lack of randomization. 

The transition from Jamming I to Jamming II takes place at (3 ~ 0.65 for all v ma x > 1- I n 
other words: When the maximum velocity is varied Free Flow, JII and coexistence phase 
(Jl+maximum current) keep constant and only the ratio between maximum current phase 
and JI changes (Figs la,b). 

V. CONCLUSIONS 

The nondeterministic Nagel-Schreckenberg model depends on the randomization, the max- 
imum velocity, and the boundary conditions: 

The buffer effect observed for the deterministic case p = and v max > 3 is strongly weakened 
with increasing randomization probability. For v max = 2 there are no buffers and therefore 
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the corresponding phase diagram is similar to the case v max = 1 for all p values. For \ max 
> 3 and p > p c (p c = 0.1172 ± 0.008 for v max = 5) the buffer effect completely vanishes 
since the development of jamming waves is no longer suppressed. As a consequence, a max- 
imum current phase occurs for p > p c and the free flow (jamming) phase can be divided 
into two regimes AI and All (BI and BII) similarly to the case of v max = 1,2. Another 
analogy to the case v max = 1 is that the free flow/jamming (free flow/maximum current and 
jamming/maximum current) transition is of first (second) order. 

There are, however, essential differences between systems with v max = 1 and v max > 1: In 
the maximum current phase the density profiles decay algebraically with an exponent 7 = 
I for v max > 2 whereas 7 = \ was obtained in the ASEP. This indicates that systems with 
v max > 1 and v max = 1 belong to different universality classes. Another difference to the 
ASEP is the existence of oscillations at the beginning of the system due to the repulsion of 
the cars. 

The comparison of systems with periodic and with open boundary conditions suggests that 
there are mainly three differences. First of all, the transition from free flow to jamming 
for systems with open boundaries is sharp and there is no maximum current phase in the 
case of periodic boundary conditions. Moreover, the dependence on the maximum velocity 
is more complex for systems with open boundary conditions due to the occurrence of the 
buffers. However, there are common features, too: Measurements of the short-range correla- 
tion function show that - as for corresponding systems with periodic boundary conditions - 
three regimes can be distinguished from each other: (a) Free Flow: cars do not hinder each 
other (b) Maximum Current and Jamming I: coexistence of freely moving and jammed cars 
(c) Jamming II: cars are jammed. 
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FIGURE CAPTIONS 

Fig la: 

Phase diagram for v max = 2 (p = 0.5, continuous line: first-order or continuous phase tran- 
sition, dotted line: border between AI/AII and Bl/Bll, broken line: border between the 
Jl and the JII regime). Although there is no symmetry along the a — (3 - line the phase 
diagram shows strong similarities to the v max = 1 case. 

Fig lb: 

Phase diagram for v m ax > 5 (p = 0.5, continuous line: first-order or continuous phase tran- 
sition, dotted line: border between AI/AII and BI/BII, broken line: border between the JI 
and the JII regime). Phase diagrams for v max = 3,4 are qualitatively the same. Due to the 
buffer effect the border between the free flow and the jamming phase shows a slight bending. 

Fig lc: 

Phase diagram for v max = 5 and p = 0, 0.1172, 0.25, and 0.5. For randomization proba- 
bilities p > p c with p c = 0.1172 ± 0.008 the buffer effect vanishes and the free flow regime 
becomes smaller with increasing randomization probability. The borders of the maximum 
current regime for a > 0.6, (3 > 0.92, and p = 0.25 (a > 0.35, (3 > 0.89, and p = 0.5) are 
represented by broken lines. The phase diagram for maximum velocities v max > 5 is very 
similar. 

Fig 2: 

Average density in the middle of the system for v ma x = 5 (p = 0.5). The first-order phase 
transition from freely moving to jammed traffic can be clearly seen. Moreover, there is a 
jump in the derivative of p(i,^) at the maximum current/free flow and the maximum cur- 
rent/jamming border what is a hint at a continuous phase transition. 
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Fig 3a: 

Current q for v max — 2, 3, • • •, 10 and (3 — 1 (p = 0.5). At a m 0.35 there is a slight 
maximum for v max > 5 which is explained in Fig 4a. 



Fig 3b: 

Current q for v ma x — 2, 3, • • •, 10 and a — 1 (p = 0.5). For v max > 5 the curves are identical. 
Fig 4a: 

Current q for p = 0, 0.125, • • •, 0.875 and (5 = 1 ( Vmax — 10). The maximum at a ~ 0.81 
for p = moves towards smaller injection rates with increasing randomization probabilities 
and for p > 0.125 a continuous phase transition is observed. 



Fig 4b: 

Current q for p = 0, 0.125, • • •, 0.875 and a = 1 (v m ax = 5). For p > 0, the buffer effect 
observed for the deterministic case vanishes and the first-order phase transition goes over 
into a continuous phase transition. 



Fig 5a: 

Density profiles for f3 — 1 (p = 0.5, v max = 5). A typical feature of density profiles in the 
free flow regime are oscillations resulting from the hindrance the cars feel at the beginning 
of the system from each other. These oscillations die out for higher system sites due to ran- 
domization. The curves decay algebraically in the maximum current regime as it is already 
known from the ASEP. 



Fig 5b: 

Density profiles for (3 = 0.7 (v max = 5, p = 0.5). The phase transition at a = 0.278 is of 
first order characterized by a linear density profile at the critical injection rate. The curve 
for a = 0.4 (a = 0.3 or a = 0.29) is typical for a density profile in the BII (BI) regime. 
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Fig 5c: 

Density profiles for a = 1 (p = 0.5, v max = 5). The results from the v ma x = 1 case are re- 
covered: In the maximum current regime (/3 > 0.89) the density profiles decay algebraically 
and in the BH-jamming regime they are described by an enhanced exponential function. 

Fig 6a: 

Correlation functions for (3 — 1 (p = 0.5, v ma x = 5). The free flow regime is characterized 
by a propagating peak with anticorrelations around it. In the coexistence regime both the 
jamming and the propagating peak can be observed. 

Fig 6b: 

Correlation functions for a — 1 (p = 0.5, v max = 5). The curves in the maximum current 
and in the Jl-jamming regime behave similarly. At f3 ~ 0.75 the propagating peak vanishes 
and the system is in the JII- Jamming ( = "Superjamming") regime. 

Fig 6c: 

Correlation functions for (3 = 1-a (p = 0.5, v max = 5). As for the injection rate a — 1, the 
propagating peak vanishes at j3 — 1-a ~ 0.75. Moreover, there is a striking similarity to 
corresponding correlation functions in the case of periodic boundary conditions. 
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